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Thermodynamics is a branch of science blessed by an unparalleled combination of generality of scope and 
formal simplicity. Based on few natural assumptions together with the four laws, it sets the boundaries 
between possible and impossible in macroscopic aggregates of matter. This triggered groundbreaking 
achievements in physics, chemistry and engineering over the last two centuries. Close analogues of those 
fundamental laws are now being established at the level of individual quantum systems, thus placing limits 
on the operation of quantum-mechanical devices. Here we study quantum absorption refrigerators, which 
are driven by heat rather than external work. We establish thermodynamic performance bounds for these 
machines and investigate their quantum origin. We also show how those bounds may be pushed beyond 
what is classically achievable, by suitably tailoring the environmental fluctuations via quantum reservoir 
engineering techniques. Such superefficient quantum-enhanced cooling realises a promising step towards 
the technological exploitation of autonomous quantum refrigerators. 

An absorption or heat-driven quantum refrigerator is a system capable of establishing a net steady-state 
transport of energy from a cold bath (c) to a hot bath (h), assisted only by the residual heat coming from an 
additional work reservoir (w) 1 " 3 . In this picture, the cold bath would play the role of the macroscopic or 
mesoscopic object to be cooled. In addition to their potential technological applications, these autonomous 
quantum- thermal devices are also appealing from the fundamental perspective, as they are naturally well suited 
for the study of thermodynamics at the level of individual open quantum systems 1,4 " 6 . 

In spite of the increasing interest that quantum absorption cooling has attracted over the last few years 5 ' 7 " 12 , the 
field is far from new. A heat-driven quantum fridge is just one specific configuration of the more general quantum 
heat pump, that can function either as a heater, a chiller or even an engine. The use of three-level solid-state masers 
as physical support for heat pumps was already discussed in the late 1950s 1314 , when spin refrigeration was also 
experimentally demonstrated 15 . The consistent quantum-thermodynamic description of these elementary three- 
level prototypes was object of further study 1 ' 4 and, just recently, alternative finite-dimensional quantum systems 
realising autonomous heat pumps have been put forward in the literature 2,3 . 

The different designs of quantum heat pumps share limitations that can be understood from the assumptions 
on their interactions with the environments. Under the familiar conditions usually met in the quantum-optical 
regime, the dissipative processes may be assumed purely Markovian 17 " 19 , which severely restricts the performance 
of any heat-driven device, and confers a distinctive spectral structure to the environmental fluctuations 16 . In 
particular, once their steady state builds up, quantum heat pumps are governed by formal analogues of the laws of 
thermodynamics and, as a consequence, their absolute efficiency ideally saturates to the corresponding Carnot 
limits s c , albeit at vanishing 'cooling power' 14 , i.e. in the reversible limit, the exchange of any finite amount of 
energy with the heat baths is performed in infinite time. 

For practical purposes, however, one needs to operate at nonvanishing power. In this case, the relevant issue to 
assess the functionality of these devices demands the optimisation of more practical figures of merit such as the 
efficiency at maximum cooling power s*. The natural question arises whether e* can approach s c arbitrarily 
closely even at finite cooling power, or if, on the contrary, it is upper bounded by some fundamental limit. The 
efficiency at maximum 'mechanical' power is extensively used to benchmark the operation of heat engines and a 
lot of effort has been devoted to establish a universal upper bound therefor 20 " 22 . Unfortunately, the general 
arguments used for engines do not provide simple bounds when the cycle is reversed into a refrigerator, and 
consequently, a different approach is needed to arrive to model-independent performance bounds. Here, we 
rigorously prove that the 'smallest' quantum absorption refrigerators, supported on ideal three-level masers, are 
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limited in their efficiency at maximum power by a fraction of e c , only 
related to the spectral properties of the environmental fluctuations at 
low frequencies. We show that this general performance bound 
applies as well to 'larger' and non-ideal designs 2 ' 3 ' 8 , as it is independ- 
ent of the details of the working material of the refrigerator. 

Achieving a good understanding of the quantum- mechanical ori- 
gin of the limitations of heat pumps can also provide key clues about 
how to surmount them. We show indeed that, by feeding an absorp- 
tion fridge with engineered thermal resources, one can push its per- 
formance bounds considerably further, allowing for classically 
impossible superefficient quantum cooling. Namely, at given fixed 
environmental temperatures, the addition of squeezing to the work 
bath leads to efficiencies above s c and, most interestingly, to a sys- 
tematic enhancement of the output harnessed power. This is 
achieved strictly within the framework of quantum thermodynamics 
and thus, in no violation of its laws 6,23 . 

Results 

Models of absorption refrigerator. As already advanced, a minimal 
model of autonomous heat pump 13 consists of a three-level system 
with each of its transitions weakly coupled to one of the three 
independent heat baths [see Fig. 1(a)]. Essentially, as the steady 
state builds up, the 'heat' collected from the cold bath is dumped 
into the hot bath with the assistance of the extra energy provided by 
the work bath, which closes the cooling cycle. Of course, the opposite 
heating cycle also takes place in the steady-state, and it is the 
imbalance between these two stationary processes which renders 
the device either a refrigerator or a heater. As we shall see below, 
refrigeration occurs as long as the frequency of the transition coupled 
to the cold bath remains below a certain threshold w c < co c ,max- 

We shall also consider the two-qubit design 3,6 of Fig. 1(b) in which 
the cold and the hot bath are each addressed through a two -level 
contact, while the work reservoir introduces dissipation in the sub- 
space {|O c )(x)|lh), |l c )®|0h)h where |0 a ) and |l a ) stand for the 
ground and excited states of contact a e {h, c}. This allows for energy 
flows between the hot and cold ends, and eventually results in net 
refrigeration. Remarkably, this fridge works within the same cooling 



window (co c 



x ) as the three-level prototype. 



Finally, we shall comment on another model 2 ' 5 ' 7 featuring three 
contact qubits connected via a three-body interaction [see Fig. 1(c)]. 
Practically- oriented issues have been recently studied in connection 
with this design, including its potential experimental realisations 910 
and the investigation of its efficiency at maximum power 8 . 



The first and second laws. When the interaction of the working 
material with the environments is sufficiently weak, one can tackle 
its effective dynamics via a quantum master equation like 

j t Q(t) = (V w + V h + V c ) Q(t). (1) 

It is an equation of motion for the reduced state of the heat pump g(t), 
where D a are the dissipation super- operators associated with each 
bath. An initial preparation p(0) = £(0)® a AfJ, factorised between 
system and baths degrees of freedom is assumed. The intrinsic 
dynamics has been eliminated by taking the interaction picture with 
respect to the free Hamiltonian of the working material H wm . 

If the dissipation is much slower than both the environmental 
fluctuations and the coherent evolution of the heat pump, the 
Born, Markov and rotating-wave approximations may be safely 
applied 16 . This leads to dissipation super-operators of the well known 
Lindblad-Gorini-Kossakovski-Sudarshan (LGKS) type 1718 , which 
are the workhorse of quantum thermodynamics 23 . We shall specif- 
ically denote them by C a in what follows: 

C a g= T coa (k(oj)qA\(oj) - X - {Al(co)A a (co),g} A (2) 

CO V / 

Here, {*, •}+ stands for an anticommutator and A a (w), for the jump 
operator associated with the decay process into channel w, which 
occurs at a rate given by the corresponding element of the spectral 
correlation tensor T^. For equilibrium reservoirs, these latter relate 
via the detailed balance condition 16 



_ -h(a/k B T a y 



(3) 



Individually, each LGKS dissipator £ a generates a completely pos- 
itive and trace preserving contractive dynamics of the working 
material, converging towards its local stationary thermal state 
T a oc exp{ — H wm /k B T a y 9 . As time goes to infinity, this contractiv- 
ity translates into (see Methods for details) 



T w Th T c 



(4) 



where T a are the equilibrium temperatures of the baths and 
Q a = tr [H wm C a g(co)] represents the energy per unit time flowing 
from bath a into the pump once in the steady state. In particular, we 
will refer to Q c as the 'cooling power' or just power. 




Figure 1 | (a) Schematic representation of the three-level heat pump. The three independent baths are taken as infinite collections of uncoupled modes 

with Hamiltonian H Ba = J^go^ The dissipative interactions in each bath-transition pair are modeled by terms like e.g. H* nt = y / y(|0)(l| + 

1 1) (0|) ®B C i with B a = g a ^ (ioc,(j. + b\ • Here, the coupling constants are proportional to y / co^, which results in flat spectral densities, and the 

parameter y controls the overall strength of the interaction 16 . The three-stroke cooling and heating processes are depicted with blue and red arrows 
respectively, (b) Diagram of a two-qubit heat pump, where we shall always assume co h > co c . Energy transport between the two contacts is mediated by 
dissipation into the work bath thanks to a coupling term of the form H 1 ^ = (| \ c ) (0/, | + |0 C ) (1/j |) ®B W . (c) A design of heat pump with a contact qubit for 
each bath. In this case, the working material is provided with an explicit three-body interaction term to allow for energy exchanges H wm =Ho+ Hi, with 
Hj= g(\ l w 0hl c ) (0 w lh0 c l + fix.). The frequencies are constrained by a>h = co c + w w 2 . 
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Since H wm is time-independent, the stationarity of the average 



energy at g(oo) implies 



Q w + Q h + Qc = 0- 



(5) 



This balance equation plays the role of a first law in quantum ther- 
modynamics as soon as one identifies the steady- state currents Q a 
with heat flows. Similarly, Eq. (4) can be regarded as a quantum- 
thermodynamic statement of Clausius theorem, i.e. the second law. 

The combination of Eqs. (4) and (5) places the ultimate ther- 
modynamic bounds on the efficiency of a heat pump in its various 
modes of operation 14 . For instance, in the chiller configuration (i.e. 
Q c > 0, Q w > 0 and Qh<0), the efficiency is defined as the ratio of 
the cooling power to the input heat (per unit time) provided by the 
work bath, i.e. e = Q c j Q w . Using Eq. (5) to eliminate Qh from Eq. (4), 
yields 



_ Qc ^ (T w -T h )T c _ 
8 ~Q W ~{T h -T c )T w -^ 



(6) 



where s c is nothing but the Carnot efficiency of a (macroscopic) heat 
driven quantum absorption fridge operating between baths at tem- 
peratures {T wy T h , T c } 24 . 

The cooling window. We now want to delimit the region within the 
space of parameters of a quantum heat pump where cooling is 
permitted by the second law. This is of course highly model- 
dependent but one can always stick to n- dimensional working 
materials with the basic three- stroke cooling mechanism of the 
three-level maser built in 13 . If the bath a couples to the heat pump 
only allowing for transitions with a gap of hco a among the eigenstates 
of H wm , and if the resonance condition w h = co c + co w holds, one 
should have 



\Qa/Qp\=CD 0L /cOp. 



This was already acknowledged as a distinctive feature of ideal three - 
level heat pumps in the seminal paper by Geusic et al. 14 , and it is easy 
to see that it remains true for the two-qubit design of Fig. 1(b). Eq. (7) 
essentially says that in a cooling cycle, every single cold excitation is 
traded for one hot excitation at the expense of consuming a single 
work excitation 7 . 

Hence, combining Eq. (4) with (7) yields 

T w 

(8) 



COu oj c 
— - + — <0 



0) c < oj c , max = 



(T w -T h )T c 
(T w -T c )T h 



€O h . 



This inequality defines the 'cooling window'. Note that the work 
temperature T w must be larger than T h (and T c ) in order to have a 
positive ft) c , max . 

In principle, the efficiency of any ideal heat pump satisfying Eq. (7) 
saturates to the Carnot bound in the reversible limit of w c — > o) Cjmax , 
i.e. when the contact transitions locally equilibrate with their corres- 
ponding baths so that the equality in Eq. (4) holds. 

On the contrary, the operation of a heat pump might become 
intrinsically irreversible if additional mechanisms of energy 
exchange were present. For instance, the three-qubit device of 
Fig. 1(c) only behaves as an ideal heat pump when the effects of 
the global interaction term Hi on the dissipative dynamics are 
entirely neglected [see caption of Fig. 1(a)] . Indeed, explicit analytical 
formulas for its steady state consistent with Eq. (7) may be written 
down in that limit 7 . However, since the contact transitions are chosen 
among the eigenstates of H 0 rather than those of H wm = H 0 + Hj, 
dissipation is always strictly 'delocalised' regardless of the interaction 
strength, and neither Eq. (7), nor the equality in Eq. (4) can be 



satisfied. Consequently, this specific design is an example of non- 
ideal heat pump, as it is unable to reach the Carnot efficiency 8 . 

Efficiency at maximum cooling power: a model-independent 
bound. In spite of its fundamental importance, the attainability of 
the Carnot efficiency in microscopic quantum heat pumps 711 is not 
the central issue for practical applications. Indeed, when operating at 
the reversible limit the power exactly vanishes. 

One would like instead to run a refrigeration cycle, carefully tun- 
ing the design parameters so that efficiency and power are maximised 
jointly, pretty much in the spirit of finite-time thermodynamics. The 
practically relevant questions would then be whether there exists a 
tight upper bound for the efficiency at maximum power g* other than 
£ c , and whether such bound is model-independent. Long-standing 
problems of this sort have been intensively studied in classical mac- 
roscopic heat engines and refrigerators 20 " 22,25 ' 26 , as well as in their 
quantum-mechanical counterparts 27 " 33 . In particular, a specific bound 
for 8* has been recently established for the model of Fig. 1(c) based 
on a numerical analysis, together with design prescriptions for its 
saturation 8 . 

In this paper, under natural assumptions on the environmental 
fluctuations, we prove analytically that the efficiency at maximum 
cooling power of any ideal fridge made up of elementary three-level 
heat pumps is tightly upper bounded by 



£*< 



d c + l 



(9) 



where d c stands for the spatial dimensionality of the cold bath. This 
applies directly to the ideal absorption refrigerators of Figs. 1(a) and 
1(b), but we verify that Eq. (9) holds as well for the non-ideal refri- 
gerator of Fig. 1(c), thus validating and generalising the bound 
obtained in 8 . Eq. (9), which is the first main result of this paper, 
establishes then a quantum-thermodynamic limitation which holds 
for all the models of quantum absorption fridges existing in current 
literature, and is manifestly independent of the details of the working 
material of the refrigerators. While full details are deferred to the 
Methods section, we sketch a proof of the bound in what follows. 

One may generically characterise the dissipation into a ^-dimen- 
sional free bosonic field a at thermal equilibrium with flat spectral 
density, by decay rates of the form 16 



r = 



y Q (D d *(\+N a (w)) w>0 
y 0 \w\ d *N a (\co\) co<0' 



(10) 



where N a (co) = [exp(h(o/kBT a ) — 1] _ 1 . Eq. (10) follows just from the 
application of the Born and Markov approximations to a general 
microscopic model of relaxation into a bosonic bath. From it, one 
can see that the cooling power of a three-level fridge [Fig. 1(a)] writes 
as 



Q c (x) =x c p(x)[aie lX — a 2 e 



-b 2 x\ 



(11) 



where x = co c /w CyinSiXy a { and b { are positive real constants, andp(*) is a 
positive function. The cooling power vanishes at both edges of the 
cooling window Q c (0) = Q c (1) = 0 and is a concave function of x. Eq. 
(11) can always be conveniently recast as 

Q c (x)=P(x) x d <(l-x)=P(x)f(x), (12) 

with P(x) again a positive function (see Eq. (41) in Methods). The 
maximum of Q c (x) is attained at some x* such that 



Q e (x,) = P'(*,)/(x») +P(x^f'(x,) = 0, 



(13) 



where the prime denotes differentiation with respect to x. Note that 
the positive and concave function f(x) is maximised precisely at x = 
d c l{d c +1). 
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Previous literature on performance bounds for heat engines has 
established that the efficiency at maximum power is maximised in 
the limit of small s c 21,34 , as it is also the case for quantum absorption 
fridges 8 . We shall therefore expand P(x) around s c ^> 0, where in fact 
P(x) ~ p(x) [see Eq. (41)], and study its analytical properties. We 
then prove the claim by reductio ad absurdum. Let us first assume 
that > d c l(d c + 1), so that the bound in Eq. (9) would not hold 
[recall Eqs. (7) and (8)] . P' (x) has at most one root in the interval x e 
[0, 1] and can be seen to be negative in the neighbourhood ofx= d c l 
(d c + 1), for vanishing e c . Being both P(x) and/(x) positive in the 
whole unit interval, the violation of the bound Eq. (9) would contra- 
dict Eq. (13). Therefore, one must have instead x* < d c /(d c + 1) and 
hence, 



co c ,. 



< 



d c 



(14) 



the bound being saturated if P'(x*) = 0. 

Now, we shall consider the two-qubit model [Fig. 1(b)]. We find 
that its steady- state heat currents can be broken up as (see Methods 
for details) 



Qc = ^c{qi+qi) 
Qh= -co h (qi+q 2 ) 



Q w = co w (qi+q 2 ), 



(15a) 



(15b) 



(15c) 



where all q t > 0. There are two contributions to the total heat fluxes, 
namely | | = w a qi and | | = oj a q 2 , which individually satisfy 
the first and second law of Eqs. (5) and (4) within the cooling window 
delimited by Eq. (8). Furthermore, | j Qp | = <J$ a / cop for i = {1, 
2} and hence, each three -level component behaves as an ideal refri- 
gerator on its own. Most importantly, the cold fluxes Q* have the 
same analytic structure of Eq. (11), so that their efficiency at max- 
imum power is limited by &^ j Sc<d c / (d c + 1). As a consequence, 

the performance of the combination of the two is also bounded by s c 
d c l{d c +1), since 



e (1 >(^)6l! ) (^)+« (2) (^)Ql 2) (^ 

qIV)+q!?(**) 



< 



d c 



(16) 



The heat currents of the three-qubit fridge [Fig. 1(c)] are much more 
involved analytically. However, extensive numerical evidence con- 
firms that Eq. (9) also applies to this case 8 [see Fig. 2(c)]. 

A numerical investigation, presented in Fig. 2, certifies that the 
bound of Eq. (9) is tight for all the three models. One can further see 
that choosing (D w /T Wf h<^l is a sufficient condition to approach it 
closely under large temperature differences T c /T/j<Cl. These are to 
be regarded as analytical design prescriptions for the practical imple- 
mentation of optimal quantum heat pumps, e.g. following recent 
proposals involving superconducting qubits or arrays of quantum 
dots 9,10 . Note as well that the limit T c /Tfc<0 also implies s* < £c^l> 
as should be expected. 

Our bound has been analytically established as a constraint on the 
performance of the fridge operation mode of the three-level maser 
and that of any ideal absorption refrigeration cycle reducible to three- 
level systems, and it has been also shown to hold for fundamentally 
different non-ideal devices. It is in order to remark, however, that Eq. 
(9) has in fact a general validity transcending specific models. Let us 
consider a totally generic model of quantum absorption fridge: We 
can always write the leading contribution to its cooling power Q c as a 
sum of terms of the form 



Q c ~hw t 



2>1 t 

CO c t - 



2>i c 



(17) 



where stands for excitation/relaxation rates of the contact tran- 
sitions at oj a . Eq. (17) just formalises the imbalance between the 
elementary three-stroke cooling and heating cycles which has to 
underly any implementation of quantum absorption cooling, and 
contains no details on the specific working material of the refriger- 
ator, nor about the spectral properties of the reservoirs to which it 
couples. Given that these three reservoirs are in thermal equilibrium, 
we can make use of the detailed balance condition of Eq. (3) and thus 
arrive to 



Qc~ha) c n&i f 



ftco c hoj w 



(18) 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



sje c ejs c ejs c 

Figure 2 | Histograms of the efficiency at maximum power e*/e c for ~ 10 5 random (a) three-level, (b) two-qubit and (c) three-qubit absorption fridges, as 
introduced in Figs. 1(a), 1(b), and 1(c), respectively. The temperatures {T a }, the hot frequency the dissipation rates, and — for the case (c)-the 
interaction strength g y were all chosen completely at random, but always respecting the conditions y^kgT^ and y<^{hco a ,g} y that guarantee the 
applicability of the master equation Eq. ( 1 ). In each case, we found the optimal pair w c } so that Q c was maximised, and computed the corresponding 
efficiency. This is equivalent to a global numerical optimisation of e* over all free parameters of the models. In each of the panels above, we show three data 

12 3 

sets that correspond to the same heat pump operating between Id, 2d and 3d baths. In agreement with Eq. (9), e* is bounded by - ec, - fie and - ec, 

respectively. It can also be seen that these bounds are tight, and that the majority of randomly sampled fridges cool very close to their corresponding 
bounds. 
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which is formally identical to Eq. (11), with the positive function 
p(x)cc Yl oc ^>^(x)/x dc ~ l . Therefore, the whole line of reasoning 
between Eqs. (1 1)— (14) is still applicable, provided that the condition 
p'(x*) < 0 is verified. This weak assumption is sufficient for Eq. (9) to 
hold, regardless of the physical support of the refrigerator and the 
specific properties of the baths. Under these premises, we conclude 
that our bound is model-independent for quantum absorption 
refrigerators. 

Let us eventually connect this result with the other known bounds 
from finite-time thermodynamics. In the seminal work by Curzon 
and Ahlborn 20 , an upper limit to the efficiency at maximum power 
was derived for a Carnot engine under the endoreversible approxi- 
mation 35 . In spite of its seemingly limited scope, such bound succeeds 
in capturing the universal behaviour of 8* at small efficiencies 21,34 , 
and recurrently appears in different models of heat engines 22,27 ' 31 ' 34 ' 36 , 
including the three-level maser in its engine configuration 34 . 
However, the obtention of similar results for refrigerators would 
require a more careful phenomenological modeling of their main 
sources of irreversibility, and can lead to highly model- dependent 
performance bounds 37 " 39 . 

In analogy with the Curzon-Ahlborn limit, our result in Eq. (9) 
accurately represents 8* at low efficiencies for all embodiments of 
quantum absorption refrigerators, which encompass fundamentally 
different models. The same limit holds as a strict upper bound to the 
efficiency at maximum power of 'classical' endoreversible absorption 
chillers 39 . It must be noted as well, that it relies on a consistent 
microscopic description of the system-baths interactions rather than 
on a phenomenological ansatz. From the physical point of view, such 
a model-independent bound can be understood from the intuitive 
notion that the maximisation of the cooling power in the steady state 
must be governed by the low-frequency dependence of the corres- 
ponding cold decay rate [see Eq. (10)]. 

Superefficiency: squeezing the second law. Up to now, we have seen 
how the standard laws of thermodynamics place fundamental 
constraints on heat pumps, even when they are made up of a 
single finite-dimensional open quantum system. This is perhaps 
not surprising taking into account that there is nothing really 
'quantum' about the operation of these devices, other than the 
discreteness of their energy spectra 8 . For instance, the ideal fridges 
of Figs. 1(a) and 1(b) operate in completely classical (diagonal) 
steady states, while the quantum coherence that builds up 
asymptotically in the non-ideal fridges of Fig. 1(c) does not seem 
to affect their performance in any crucial way. Even if steady- state 
bipartite entanglement may exist in this case, it usually appears only 
in the regime of very low efficiencies and vanishing cooling power 11 . 
Likewise, other types of quantum correlations, though widely 
present, do not have any influence on optimal cooling 8 . 

Is it then possible at all for quantum heat pumps to operate past the 
'classical' limits? Ideally, one would like to devise tricks to push their 
performance bounds further, of course remaining always within the 
standard framework of quantum thermodynamics and not advoc- 
ating any violation of its laws. We shall devote this section to illus- 
trate how one can indeed go beyond Eqs. (6) and (9), by exploiting 
non-equilibrium environmental features that can be mimicked with 
suitable quantum reservoir engineering techniques. 

The whole idea would consist in initialising the work reservoir in a 
squeezed-thermal state 

xy = S{H)Xl${ti). (19) 

Such a state results from the action of the unitary squeezing oper- 
ator 40 

S(f)=IIexp|ir^-if(fet) 2 J (20) 



on a thermal preparation X^, where |^,^| stand for creation and 

annihilation operators on mode \i of the work reservoir. The para- 
meter £ = r exp iO is in general a complex number, although for our 
purpose we can set the phase 0 to zero and consider a real squeezing 
parameter £ = reU. The quantumness of a squeezed state resides in 
the asymmetry of the variances of its field quadratures X+ =b f/L ±P fl , 
as opposed to a (classical) coherent state, for which AX+ = AX^_ . 

It is well known that squeezed states, even in absence of entan- 
glement, allow for quantum enhancement in several applications 
of information theory 41 , including quantum cryptography 42 , and 
most notably precision measurements 43 and quantum metrology 
below shot noise 44 , with applications e.g. for gravitational wave 
interferometry 45 47 . 

The key property to be exploited here is the non-stationarity of 
squeezed preparations. This results in a periodic time modulation of 
the correlation functions of the work reservoir, which is somewhat 
equivalent to the action of an external driving on the heat pump. In 
this sense, the work reservoir may now play an active role in the 
cooling process. As we shall now see, this is achieved without com- 
promising the thermodynamic consistency of the whole setting. 

In order to account for the squeezing of the work reservoir, the 
quantum master equation Eq. (1) must be replaced by 

j t e(t) = (& r w +£h+C c )e(t), (2i) 

where the modified work dissipator j£f r w = C r w + AC r w can be cast in 
the standard LGKS form 16 (see Methods for details). 

Clearly, if the work heat current is suitably redefined as 
Q w = tr[H wm j5fJ v ^(oo)], the asymptotic stationarity of energy 
implies again the first law as stated in Eq. (5). 

However, the dissipator Jjf r w will not generally bring the corres- 
ponding contact transition to thermal equilibrium at temperature 
T w , but to some other steady state o w (r). Still, the LGKS form of 
the overall generator of the dynamical map in Eq. (21) guarantees its 
full contractivity 19 , so that the second law generalises to 

-trfe(«D)log^(r)] + ^ + ^<0. (22) 

It is convenient to fit o w {r) by a thermal state of the form 
o w (r) oc exp{ —H wm / k B T^ { (r)}, with a suitable squeezing- depend- 
ent work temperature T^ ff (r). Since the contact with the work res- 
ervoir only involves two levels, this is always possible. Eq. (22) maybe 
thus rewritten in the more familiar way 

Q w ,Qh,Qc <0 (23) 

For any r > 0, the effective temperature T^ ff (r) exceeds T w , and it 
diverges for r — » °° Therefore, the generalised squeezed-dependent 
Carnot efficiency will always exceed its classical value 

£cW ^( 1 -5?k) >£c(0) - (M| 

Note that, as should be expected, in the limit of r — » o°, Eq. (24) 
saturates to the maximum efficiency of a power- driven refrigerator 
(i.e. a reversed heat engine) £c(°°) = T c /(T h — T c ) 14 . 

Hence, we have demonstrated how, for a given set of thermal 
resources {T c , T h , T w }> classically forbidden superefficient quantum 
absorption cooling may be realised by just squeezing the work res- 
ervoir. This is illustrated in Fig. 3, with a set of performance char- 
acteristics corresponding to different environmental squeezing 
parameters in the three-level maser heat pump of Fig. 1(a). 
Perhaps even more remarkable than the seemingly counter-intuitive 
possibility of cooling above the classical efficiency threshold, is the 
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Figure 3 | Parametric plot of power vs. efficiency of a three-level maser 
operating between three-dimensional bosonic reservoirs at temperatures 
T w = 170, T h = 80 and T c = 30 (h = k B = 1), for different values of the 
squeezing parameter of the work reservoir r = {0, 0.5, 1, 1.5, 2}. Both Q c , e* 
and the maximum e increase with r. In each curve, the hot frequency was 
fixed at co h = 50, while co c ranges from 0 up to co C)max (r). The frontier 
between the white and the shaded region corresponds to the performance 
characteristic of an ideal reversed heat engine operating between the same 
hot and a cold baths. The general behaviour illustrated here is totally 
independent of the specific values of and the set of temperatures { T a }. 

fact that the efficiency at maximum power and the cooling power 
itself systematically increase jointly with the squeezing parameter. 
This is a striking consequence of the non- equilibrium environmental 
fluctuations. 

In the macroscopic domain, the practical drawbacks of an absorp- 
tion chiller 48 , as compared to an equivalent work- driven refrigeration 
cycle, are typically its much lower efficiency and output power. 
Nevertheless, as we can see in Fig. 3, a quantum-enhanced heat- 
driven fridge does not only outperform its purely thermal counter- 
part, but is also capable of approaching very closely the performance 
of a Carnot fridge, when provided with a finite amount of envir- 
onmental squeezing. Thus, as opposed to the classical case, quantum 
absorption cooling may compete on an equal footing with conven- 
tional refrigeration, when combined with reservoir engineering tech- 
niques 49 " 51 to allow for the exploitation of non-equilibrium squeezed 
environmental fluctuations. This observation could be a promising 
first step towards a number of applications of heat-driven refrigera- 
tion to a new generation of quantum and nanoscale technologies. 

Discussion 

Let us start by summarising the two main results of this paper. On the 
one hand, we proved that the efficiency at maximum power of all the 
known models of quantum absorption refrigerator (both ideal and 
non-ideal) is tightly upper bounded by a fraction of the Carnot 
efficiency £ c , which is independent of the details of the device and 
only relates to the spectral properties of the thermal fluctuations of 
the environment. On the other hand, we showed most remarkably 
how by squeezing the heat source that drives an absorption cooling 
cycle, one may boost its performance to the extent of making it 
comparable to a conventional power- driven cooling device. 

A key path to the optimisation of autonomous quantum heat 
pumps is thus found to reside in applying suitable reservoir 



engineering techniques 49 " 51 rather than in exploiting any resource 
intrinsic to their quantum mechanical working materials 811 . 
Indeed, the importance of reservoir manipulation has been very 
recently acknowledged in the context of quantum thermodyn- 
amics: It has been recently speculated that super efficient opera- 
tion of quantum heat engines, in apparent violation of the second 
law, may be achieved e.g. by reservoir squeezing 52 or using more 
general types of non -equilibrium reservoirs 53,54 , and by connect- 
ing the working material to an adiabatically isolated ancilla 55 . 
However, to our knowledge, this paper represents the first 
demonstration of a systematic enhancement in the performance 
of a quantum-mechanical thermal device. Furthermore, this 
unconditional enhancement is, in principle, achievable without 
manipulating the working material or the given hot and cold 
baths: It is only the arguably controllable heat source that needs 
to be tailored. 

At this point, one could wonder if generating squeezing in the 
work reservoir is really worth the effort when one could raise instead 
its equilibrium temperature to obtain effectively the same results. 
Back to Fig. 3, we see that a squeezing parameter of r = 1.5 
(—13 dB, which is currently at reach 45 ), would already take the heat 
pump close to its best equivalent power- driven counterpart. In order 
to achieve a comparable amplification without squeezing, one should 
increase the work temperature by at least one order of magnitude, 
which might just not be possible in engineered or natural environ- 
ments. Thus, it seems reasonable that in many concrete situations, 
reservoir squeezing could be indeed the best and most natural way to 
boost the performance of a heat- driven fridge. However, supporting 
the claim that a quantum- enhanced absorption refrigerator can 
really compete with a heat pump driven by work is a more delicate 
issue: To mimic a non -stationary state for the work bath 49 " 51 , one has 
pay an extra cost that should be added to the overall efficiency bal- 
ance for a fair comparison. Again, depending of the specifics of the 
implementation, it may well be that a quantum- enhanced quantum 
absorption chiller came to cool cheaper than the corresponding 
quantum 'compression' cycle. Furthering this issue demands a study 
on its own. 

We also wish to make clear the realm of validity of our analysis. In 
all of the above, we have taken the fulfilment of the first and second 
laws as a guarantee of thermodynamic consistency, but we have 
neither discussed the third law, nor probed the neighbourhood of 
the absolute zero 6 . Since we approach the open system dynamics 
from a LGKS-type master equation, the Markov approximation 
should better hold, which in turn implies that the thermal fluctua- 
tions must be sufficiently fast as compared with the dissipation time 
scales. Thus, to be always on the safe side, we must limit our envir- 
onmental temperatures from below. 

When it comes to the experimental realisation of quantum- 
enhanced absorption technologies it is worth mentioning the pio- 
neering refrigeration experiments by Geusic et al. 15 using the energy 
level structure of Cr 3+ ions in a ruby crystal as support for a three- 
level maser. As already mentioned, detailed proposals exist as well for 
the non-ideal three-qubit model of Fig. 1(c) using superconducting 
qubits 9 , and quantum dots 10 . The squeezing can be engineered by 
coherently driving the work transition, potentially involving auxili- 
ary levels and vacua; feasible schemes have been proposed e.g. invol- 
ving trapped ions or Rydberg atoms 49 " 51 . 

In conclusion, we have demonstrated the possibility of genuine 
quantum-enhanced absorption refrigeration beyond the fun- 
damental bounds imposed by classical thermal environments. The 
application of reservoir engineering techniques to autonomous 
quantum heat pumps might render them practically useful and com- 
petitive for many applications of quantum technologies, in primis 
quantum cooling. Moreover, the distinctive simplicity of these 
devices is ideal to get a clean glimpse of how thermodynamics looks 
like beyond the standard scenario. 
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Methods 

Quantum master equation. We will now write down explicitly the equation of 
motion for the working material of a quantum heat pump, introduced in Eqs. (1) and 
(21). The total Hamiltonian is generically (see caption of Fig. 1) 



H — H wm + H Bx + 



(25) 



where the dissipative system-reservoir interactions H int may be taken as 

H^x(|0 a >(l a | + |l a >(0 a |)®B a s6- a ®B B . (26) 

Here, the ground and excited states of the two-level contact port with bath 
a are denoted by |0 a ) and |l a > respectively. Let us first assume that all three 
reservoirs are prepared in thermal equilibrium. Under the further assumptions of 
vanishing initial correlations between system and environments, weak system- 
environment interaction and separation of time scales of free and dissipative 
dynamics (Born, Markov and rotating-wave approximations), one arrives to the 
well known LGKS-type quantum master equation 16 with dissipators as those of 
Eq. (2). 

The elements of the spectral correlation tensor, given explicitly in Eq. (10), 
follow from the power spectrum of the environmental correlations 

r ro!a = 2Re| ds e im (^Bl(t)B a (t-s)y On the other hand, A x (co) result from the 

decomposition of a a as eigenoperators of H wm 16 , and the discrete index co labels all the 
open decay channels, i.e. all the energy differences that correspond to non-vanishing 
jump operators. 

In the case of an ideal fridge, for which |0 a ) and | l a > are picked among the eigen- 
states of H wm , Eq. (1) rewrites as 



dg 
~dt~~ 



Y^ T *><°\°* QK~2^ d oi >£}+ + 

a cogQ ^ 



e -*°>/k.T. i +~- 



(27) 



where <7+ = | l a ) (0 a | and o~ = (d + and the inner summation runs over the 
frequencies Q = {co w , co h , co c }. Eq. (27) accounts for the reduced dynamics of the 
three-level and two-qubit heat pumps of Figs. 1(a) and 1(b). The case of the non-ideal 
three-qubit model of Fig. 1(c) is more involved. All the details may be found 
elsewhere 8 . 

Let us now relax the assumption of equilibrium environments to allow for 
squeezing in the work reservoir. The key difference is that the rates T a>(0 become 
explicitly time-dependent as a consequence of the non-stationarity of squeezed pre- 
parations. After performing the rotating-wave approximation, new terms appear in 
the work dissipator 16 , that becomes 

J&T W = C r w + AC r w = £ r U UwH QAl H - i { Al(cD)A w (co),g} + ) 

0) \ / 

(28) 

-A w (co) 2 q-qAI(co) 2 ), 
where the squeezing-dependent coefficients Y r w m and A r w 0J are given by 

' y 0 co dw (N r w (co) + l) co>0 



YoW dw K(\cn\) w<0 



(29a) 



^tr[(^)log^] + -S(^)>0. 



(32) 



where S(g) = —tr[g log g] stands for the von Neumann entropy. We may replace g in 
Eq. (32) with the steady state of the full dynamics g( go ), and the local steady states g° a , 

with equilibrium states e~^™A BT «. Note that indeed C a e~ Hv,m l k]iTa =0 19 . We thus 
obtain 



^ tr[H wm Ag(oo)] >Q ^ Q w | Q h | g £<Q 



k B T a 



T w T h T c 



(33) 



that is, we recover the second law as stated in Eq. (4). When the work reservoir is 
prepared in a squeezed thermal state, the derivation remains exactly the same, only 
replacing T w with the squeezing-dependent effective temperature 52 



hco w 



tanh r-\-expHco w /k B T w 



> T 



(34) 



k B \og 



1+tanh rexp Hco w /k B T w 



Tjf (r) is such that ^ e - H ™/W ff ( r ) =0. This leads to the modified second law of 
Eq. (22). 

Efficiency at maximum power of a three-level fridge. We will devote this section to 
prove the ultimate bound s* < d c /(d c + l)s c on the efficiency at maximum power for 
the three-level prototype of absorption refrigerator of Fig. 1(a). From Eq. (27), the 
steady state of the three level system may be readily found to be 



011(OO): 



£22(00) = 



£33(00) 



v +r_ (Uc r COw 



(35a) 



(35b) 



(35c) 



Here, the ground, first and second excited states of H wm are labeled 1 1), |2> and |3> 
respectively and all coherences asymptotically vanish. T Wa is a shorthand for r a (0fl 
and the denominator A is given by 



a = r^r^ + r mw + r_ m r Q}w + r^r^ + r^r* 
+r_ ffl r mh +r_ ffl ,r ffl +r_ co ,r_ f0 +r_ f0 r_ fflw . 



(35d) 



Combining the excitation and relaxation rates r_ ffl(: , of the cold transition with 
its steady-state populations g n (°°), g 22 (°°) one obtains a stationary heat current Q c 
given by 



Q c = co c ( r _ 0Jc g l ! ( 00 ) - T 0Jc g 22 ( 00 ) ) 



(36) 



which, using the detailed balance relations, becomes (from now on we will set H = k B 
= 1) 



Qc=co c r Wh r (0w r t 



e -co c /T ce -co w /T w _ e ~coh/T h 



(37) 



- l -yo\co\ d »M r w (\co\), 



(29b) 



withiV^(a»)=iV w (ft»)cosh2r+ sinhr 2 and M r w = smh2r(2N w (co) + l)/2. A new set 
of jump operators can always be found so that the work dissipator in Eq. (28) takes the 
standard LGKS form of Eq. (2), which in turn, guarantees that the effective dynamics 
of the working material remains that of a dynamical semigroup (i.e. completely 
positive and trace-preserving) in spite of the squeezing 17 . 

The second law at steady state. It is known that the entropy production of a 
dynamical semigroup is a strictly positive quantity 19 . It is defined as 

j2m)]^-j t s(m\\~e°)>o, (30) 

where S(£||£°) =tr [£(log g — log g 0 )] is the quantum relative entropy and £° is a 
steady state of the dissipative dynamics. Each of the dissipators appearing in Eq. (1) 
individually satisfies the inequality 

tr[(£ a £)log$] -tr[(£ a £)log£]>0, (31) 

where g° a is now stationary to the local dissipator £ a alone, i.e. C a Q Q a = 0. Now, 
summing Eq. (31) over a yields 



Given the constraint co^ = co c + co w , we may fix the values of {co w , T w , T^, T c } and take 
x = oj c / (co w s c ) as the only independent variable in Eq. (37). Of course, one would 
obtain the same results by fixing co h instead of co w . This brings us back to Eq. (11) 



Q c (x) =x dc p(x) (a x e 



(38) 



where a\ =e~ 0)w//Tw , a 2 = e~ 0Jw ^ Th , b x =co w s c /T c , b 2 = co w 8 C /Th andp(x) is positive. 
As already mentioned, the cold heat current Q c (x) is concave and positive in the unit 
interval, vanishing at the boundaries Q c (0) = Q c (l) =0. 

Let us take a closer look to g(z)=z dc [a\e~ hlZ — a 2 e~ hlZ ). Its analytic continuation 
into the complex plane has, in addition to its zeroes on the real axis, the following 
complex roots 

2%i T h T w 



Zn = l+ 'W 1 '^^ n ^\{0}. 

co w {T w -T h ) 



(39) 



By applying the Hadamard factorisation theorem 56 , one may conveniently rewrite 



sinh 



g{z)=z d < ae- cz {l-z)- 



T w -T h 



(1 —z) sinh 



co w T w 



2 T w T h 



(40) 
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with a,ceR* 
becomes 



, In particular, c = co w s c (T c + T h )/2T c T h . Back into the real axis, Q c (x) 



Qc(x) 



ae cx p(x)smh e(l— x) 
(1 — x)sinh e 



Hi- 



(41) 



that is, we recover Eq. (12) by identifying the first factor in the r.h.s with the positive 
function P(x). The dimensionless constant e is defined as e= (T w — Th)co w /2T w Th. 
From there, the performance bound e* < d c l(d c + 1) follows without difficulties [see 
Eqs.(12)-(14)]. Note that, the fact that P'(x) has at most one zero in the unit interval 
follows from the concavity of Q c (x). 

From Eq. (41), one sees that taking s c ^> 0 yields P(x) ~ p(x), while p(x) itself tends 
to the constant value Q c (x*): That is, Q c (x) / Q c (**) ~x dc (l—x) =f(x). To probe 
this limit, one can just expand P(x) around e c ^ 0 and take its derivative with respect 
to x at x = d c /(d c + 1). This yields 



p, wL= 



= cosh e— cosh 



e(l-d c ) 



d c + l 



2d c e sinh e 
' {d c + lf 



+ 0 «*« 



(42) 



which is a strictly negative function of e. Consequently, P' (x) < 0 at x = dj(d c + 1 ) for 
s c — » 0. Indeed, it can be seen that P'(x) < 0 e (0, 1) in this limit. 

Three-level breakup of the two-qubit fridge. We shall finally show how the two- 
qubit fridge of Fig. 1(b) can be broken up into two coupled three-level masers. As a 
consequence, following the steps of the preceding section, one can prove that the 
performance bound of Eq. (9) also applies to the two-qubit refrigerator. 

Let us denote the eigenstates of the working material H wm by { |0 h 0 c ), |0 h l c ), | lhO c >, 
|lhl c )}> with energies {0, co c , co h , co h + co c }. The cold bath couples locally to the cold 
qubit [see Fig. 1(b)] and thus, it only drives the transitions |0 h 0 c ) <-> |0 h l c ) and | l h 0 c ) 
<-> | l h l c >. On the other hand, the hot bath is connected to |0 h 0 c > <-> | l h 0 c > and |0 h l c > 
<-» |lhl c )> while the work reservoir operates in the subspace |0 h l c ) <-> |lhO c } [see 
Fig. 4]. 

Intuitively, one can already see that the two-qubit fridge is comprised of two 
elementary three-level masers sharing the work transition. To make this statement 
more precise, let us write down explicitly the stationary heat currents 

Q c = fflc(<O h l c |£ c §(cx))|O h l c ) + (lhlc|/: c §(a))|lhlc)) 
Qfc = Q)h«lhO c |Ae(oo)|lhO c > + (lhl c |Aff(oo)|lhl c » 
Q w = co w (l h 0 c \C w g(co)\l h 0 c ), 

where co w = co h — co c (co h > co c ). From here, we can recover Eqs. (15) by just defining 

qi = (o h i c |£ c ^(cx))|o h i c > = - <ihO c |Ae(oo)|i h o c ) 

g 2 = (lhl c |£ c ^(cx))|l h l c ) = -(l h l c |£^(cx))|lhlc> 
q\ +qi = <lh0 c |/: w ^(oo)|l h 0 c >. 



As already pointed out, within the cooling window co c 



, the heat currents of 



each of the two three-level components Q a behave as those of two ideal fridges, 
individually satisfying the first and the second laws of Eqs. (5) and (4). Closed 

formulas for and Q,f\ analog 
state solution of Eq. (27), namely 



formulas for and Q,f\ analogous to Eq. (37), maybe obtained from the steady- 



where 



-7 l t r 1 

-^01 A 0) c , A _ 



(44a) 



llfclc) 

|i fc o c > 

|0fcl c > 

\0 h 0 c ) 



I 



I 



I 



Figure 4 | Level diagram of the two-qubit fridge of Fig. 1(b) with the 
transitions coupled to the work, hot and cold reservoirs depicted as red, 
orange and blue arrows respectively. 



and 



f 1 

(Oh 


— 7 l T 


r 1 i 

1 -CO h — 1 -<*>h 


(44b) 






+a 1 r (0;i r_ co j 


(44c) 






Jw +air_ aj/i r COc ) 


(44d) 




=2 0 1 ir_ (De 




(45a) 




= Z\()F -CQ h 




(45b) 


f 2 


= Z\o(y -co 


,+a 2 r_ Q , fc r a)c ) 


(45c) 


f 2 

— CO 


„ = ^01 (^« H 


+a 2 r a;ft r_ (0 j. 


(45d) 



The denominator A z has the same form of Eq. (35d), only replacing T 0h with f l 0J 
The remaining constants in Eqs. (44) and (45) are 



^r^+r^)- 1 , a 2 = (r 



r^+r^+r. 



7 1 - 

' ^10 ~ 

- 7 2 

5 1 ( 



+r fflc +r_ 



+r_ COc +r £0 



It is clear that both and can be cast in the generic form of Eq. (38), so that 
their associated efficiencies at maximum power must be bounded by 

^<d c /{d c + l). 
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